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We consider two-dimensional non-centrosymmetric superconductors, in which the order parameter 
is a mixture of s-wave and p-wave parts, in the presence of an externally induced Zeeman splitting. 
We derive the conditions under which the system is in a non-Abelian phase. By considering the 
non-degenerate zero-energy Majorana solutions of the Bogoliubov-de Gennes (BdG) equations for a 
vortex and by constructing a topological invariant, we show that the condition for the non-Abelian 
phase to exist is completely independent of the triplet pairing amplitude. The existence condition 
for the non-Abelian phase derived from the real space solutions of the BdG equations involves the 
Pfaffian of the BdG Hamiltonian at k — 0, which is completely insensitive to the magnitude of the p- 
wave component of the order parameter. We arrive at the same conclusion by using the appropriate 
topological invariant for this case. This is in striking contrast to the analogous condition for the 
time-reversal invariant topological phases, in which the amplitude of the p-wave component must be 
larger than the amplitude of the s-wave piece of the order parameter. As a by-product, we establish 
the intrinsic connection between the Pfaffian of the BdG Hamiltonian at k — (which arises at the 
BdG approach) and the relevant Z topological invariant. 



PACS numbers: 03.67.Lx, 71.10.Pm, 74.45.+C 

In quantum mechanics, statistics is defined as the 
transformation rule for a many-particle wave function 
under a pair-wise interchange of the particle coordinates. 
Bosons and fermions are the simplest of the quantum par- 
ticles in that the corresponding many-body wave func- 
tions undergo a trivial transformation - multiplication 
by 1 or -1, respectively, - under an interchange of the 
quantum numbers of any two particles. In the special 
case of (2 + 1) dimensions, in which simple permutation 
of the particle coordinates and actual exchange opera- 
tions can be inequivalent, the statistics of particles can be 
more complex. [TH3] Anyons are those quantum particles 
for which the two-dimensional many body wave-function, 
under pair-wise exchange of the particle coordinates, re- 
ceives a phase factor e 10 , where the statistical angle 9 
can take any value between (bosons) and ir (fermions) . 
An example of anyons are the quasiparticle excitations 
of the celebrated v = | fractional quantum Hall state in 
which the quasiparticles have anyonic statistics with the 
statistical angle 9 = | . [I] 

In even more complex two-dimensional systems, if the 
many-body ground state wave function happens to be 
a linear combination of states from a degenerate sub- 
space, a pair-wise exchange of the particles can unitarily 
rotate the wave function of the ground state in this sub- 
space. Consequently, the statistics is non-Abelian, |6] 
and the corresponding system is a non-Abelian quantum 
system. Non-Abelian quantum systems in the so-called 
Ising topological class, [6] are characterized by quasipar- 
ticle excitations called Majorana fermions. In supercon- 
ducting systems, which are our focus in this paper, Ma- 



jorana fermions arise as non-degenerate, spatially local- 
ized, zero-energy quasiparticle excitations bound to de- 
fects (such as vortices and boundaries) of the supercon- 
ducting order parameter. Because of the non-degeneracy, 
the Majorana fermion solutions of the excitation spec- 
trum are topologically protected, i.e., any local pertur- 
bation to the BdG Hamiltonian near an order param- 
eter defect cannot move a Majorana fermion solution 
away from zero energy. The second quantized opera- 
tors, 7^, corresponding to these zero energy excitations 
are self-hermitian, 7I = 73, which is in striking contrast 
to ordinary fermionic (or bosonic) operators for which 
Cj 7^ c\. The Majorana fermions, which are actually more 
like half-fcrmions, were envisioned [7] by E. Majorana in 
1935 as fundamental constituents of nature {e.g. neu- 
trinos are thought to be Majorana, rather than Dirac, 
fermions). Majorana modes are intriguing [5] because 
each Majorana particle is its own anti-particle unlike or- 
dinary fermions where electrons and positrons (or holes) 
are distinct. 

There has been a growing interest in realizing Majo- 
rana fermions in the laboratory for topological quantum 
information processing purposes. 6J Majorana fermions 
can be used to build topologically protected quantum 
memory applications. Even more importantly, they can 
be used in topological quantum computation (TQC) 
along with supplementary unprotected quantum gates re- 
quiring only small amounts of error corrections. [9 TQC, 
in contrast to ordinary quantum computation, would not 
require any quantum error correction since the Majorana 
excitations are immune to local noise by virtue of their 
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non-local topological nature. [SJ [TU] One of the first can- 
didate systems supporting Majorana fermion excitations 
that have been proposed are the non-Abelian v = 5/2 
fractional quantum Hall state [11] and chiral p-wave su- 
perconductors [T3] or supcrfluids. Recently, Fu and Kane 
[13] proposed that the surface of a strong topological 
insulator (STI) in proximity to an s-wave superconduc- 
tor can support a single non-degenerate Majorana mode 
bound to the vortex cores. Later, Sau et al. [H] sug- 
gested that the Fu and Kane set up can be significantly 
simplified by replacing the STI with an ordinary semi- 
conductor with strong spin-orbit coupling. It was shown 
that in the presence of a proximity induced pure s-wave 
pairing potential A s and a Zeeman splitting, the sys- 
tem under certain conditions [141 115] can support a non- 
degenerate Majorana fermion in the vortex core. Fol- 
lowing this, it was quickly realized [16] that the one- 
dimensional version of the same set up, a semiconduct- 
ing Majorana nanowire with zero-energy states at the 
two ends, would likely be an easier system to explore the 
physics of Majorana fermions, since the gap to the higher- 
energy non-topological excitations (so-called mini-gap) is 
of order A s (there are no other sub-gap states other than 
the Majorana states). Further, the Zeeman splitting in 
this case can be introduced by a magnetic field parallel 
to the superconductor, |17U18j making a proximate mag- 
netic insulator |14j unnecessary. A comprehensive discus- 
sion of the spin-orbit coupled semiconductor with prox- 
imity induced s-wave superconductivity has been given 
in Ref. [15]. 

Even though spin-orbit coupled systems with proxim- 
ity induced s-wave superconductivity have some impor- 
tant advantages, such as mini-gap ~ A s much larger 

A 2 

than the usual mini-gap ~ (A< is the p-wave gap 
and Ep is the Fermi energy) in chiral-p-wave supercon- 
ductors, a crucial new requirement is good supercon- 
ducting proximity effect itself. Therefore, it is impor- 
tant to explore the possibility of Majorana fermions in 
other spin-orbit coupled systems which are intrinsically 
superconducting, obviating the need for the proximity 
effect. In this paper, we explore this possibility in non- 
centrosymmetric superconductors (NCS) where the su- 
perconductivity is inherent to the material. Gor'kov and 
Rashba have shown [20] that strong spin-orbit interaction 
near a doped surface of a three-dimensional BCS-type su- 
perconductor lifts the two-fold spin degeneracy and the 
lack of inversion symmetry at the surface results in a 
two-dimensional pair wavefunction which is a mixture of 
both s-wave and p-wave parts. Penetration depth stud- 
ies on NCS like Li2Pd 3 B and Li2Pt 3 B with purely s-wave 
interactions, have shown that the spin-orbit interaction 
and lack of inversion symmetry result in s-wave and p- 
wave mixing. |21j The superconducting transition tem- 
perature (T c ) for Li 2 Pd 3 B is 6.7 K and that for Li 2 Pt 3 B 
is 2.43K. |21j Similar admixture of s-wave and p-wave 
Cooper pair wavefunctions have been also been predicted 
in heavy fermion compounds like CePt 3 Si with broken in- 
version symmetry. [22) A NCS with an externally applied 



Zeeman splitting is thus another candidate for realizing 
non-Abelian phases and Majorana bound states (MBS) 
in vortex cores. 



However, due to the mixing of s-wave and p-wave 
pairing in NCS, it is apriori not obvious that the con- 
dition for the existence of Majorana modes as derived 
in Ref. Q3] would hold even for NCS. Note that this 
condition, namely, Vf > A 2 + /i 2 where V Z ,A S and fi 
are the Zeeman splitting, s-wave pair potential, and the 
chemical potential, respectively, is derived in Ref. [14] 
for the case when the p-wave pair potential is absent. 
In that case, since the superconducting pair potential is 
taken to be proximity-induced by an adjacent s-wave su- 
perconductor, 13] the question of a p-wave part of the 
pair potential does not arise. In this paper we will con- 
sider the spin-orbit coupled systems when both s-wave 
and p-wave pairing potentials are simultaneously present. 
While in some of the real systems the s-wave part of the 
pair potential has been found to be dominant over the 
p-wave component, in some other systems the p-wave 
part is dominant [23[ . By calculating the appropriate 
topological invariant and explicitly solving for the MBS 
using BdG equations, we shall show that exactly the 
same condition as stated above still holds for the exis- 
tence of the non-Abelian state even in the presence of 
an explicit p-wave pair potential. This is important be- 
cause it means that, irrespective of the relative values 
of the s-wave and the p-wave parts of the order param- 
eter, a non-centrosymmetric superconductor can always 
be brought into a non-Abelian phase by the application 
of a sufficient Zeeman splitting. In contrast, existence 
of Majorana fermions in non-centrosymmetric supercon- 
ductors was found earlier [24, 25j only for the case when 
|A t /A s | > 1, where A t and A s are the p-wave and the 
s-wave parts of the total pairing potential. 



The paper is organized as follows: In section I, we shall 
write down the BdG Hamiltonian for a 2D NCS with 
Zeeman splitting. We shall then give an explicit formula 
for the topological invariant and derive the condition for 
which the system is in a topologically non-trivial phase. 
In section II, we shall consider the non-degenerate Ma- 
jorana bound states in vortex cores of such a NCS with 
Zeeman splitting. By writing down an effective ID BdG 
equation and considering its solution out side the vortex 
core, we shall arrive at the condition for the MBS to ex- 
ist. In section III, we shall consider the edge states in 
such a system and show that a non-degenerate MBS in 
vortex cores necessarily leads to gapless edge states. In 
section IV, we shall demonstrate the connection between 
the invariant in section I and the Pfaffian at zero mo- 
mentum of the BdG Hamiltonian considered in section 
II. We shall summarize our findings in section V. 
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I. CALCULATION OF TOPOLOGICAL INDEX 
FOR NCS WITH ZEEMAN SPLITTING 

In this section we give an explicit formula for the topo- 
logical invariant for non-centrosymmetric superconduc- 



tor (NCS) in 2D and in the presence of a Zeeman split- 
ting. The Hamiltonian for a two-dimensional NCS with 
both s-wave and p-wave pairing amplitudes and a Zee- 
man splitting applied in the z direction has the following 
form 



£ fe <7 + V z a z + a(ar x k).z A s a + A t (cr x k).l 
A s o- + A t (<T x k).z £ k <r + V z a z - a(cr x k).£ 



Ck 

J 



(1) 



Here £ fc = ri(k 2 x + k 2 y ) - fi, 77 = h 2 /2m, (c£, c_ fe ) = 

( c l t c k l c -k,l ~~ c -fc,t)' Vz is a perpendicular Zeeman 
splitting which can be externally induced, a is the spin- 
orbit coupling, and A s and A t are the s-wave and the 
p-wave components of the pair potential, respectively. 
We note that the BdG Hamiltonian without time rever- 
sal and spin rotation symmetry belongs to class D [26] 
and is characterized in 2D by a topological invariant C\ , 
which is the first Chern number of the £7(1) bundle de- 
scribing the many-body wavefunction. The topological 
index C\ written in terms of the quasiparticle Green's 
function G(k) has the following form: 

Ci = ^2 J #lsdu(?t[Gd km G- 1 Gd ky G- x Gd u G-' L ] 

- Tr [Gdk y G~ 1 Gdk x G~ 1 Gd^G^ 1 ]) (2) 

Volovik and Yakovenko have derived the same expression 
for the topological invariant in He 3 -A [27]. Writing the 
inverse Green's function as G -1 = H — iui, and using 
Eq. (JlJ we numerically compute Gi using the formula ([2| 
for a range of parameters. We find that the system is 
topologically non-trivial (C\ = 1 with error bars ±2 x 
10~ 8 ) as long as the following condition is satisfied: 

V 2 > v 2 + A 2 S (3) 

Note that this is the same constraint as obtained in 
Ref. [14] for the case when the s-wave pairing potential 
is proximity induced and, therefore, a triplet pair poten- 
tial is absent. We thus find that the condition for the 
existence of non-Abelian phase remains unchanged even 
in the presence of triplet pairing and is thus completely 
independent of the ratio of triplet and singlet pairing am- 
plitudes. 

This result is a major improvement over Ref. 24 where 
only the case |A t /A s | > 1 has been considered. In that 
work it has been claimed that, in order for the non- 
Abelian phase to be stable, the single particle energy 
gap should not close. A sufficient condition for the single 
particle gap to not close is |A t /A s | > 1. This condition 
for the existence of the non-Abelian phase is required if 
the Chern number is calculated (as is done in Ref. |24] ) 
by decoupling the spins in the limit a — ¥ 0, A s — > and 
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FIG. 1: Plot shows how the energy gap scales with the ratio 
A t / A s for V z = 6.0, fi = 4.0, a — 1.0, and different values 
of A s (all energy scales are in units of 77 = h 2 /2m). The 
discontinuity in the gap for A s = 3.0, A s = 2.0, and A s = 1.5 
are due to the global minimum of the gap function switching 
from k ^ to k — (see inset of Fig. [2j|. Eventually, for very 
large values of A s , the global minimum always sits at k = 
and there is no gap closing at any value of At/A s . 



mapping the system to two independent copies of spin- 
less p x + ip v superconductors. One can then calculate 
the TKNN numbers [55] for the two spins independently 
provided the gap does not close during the mapping, and 
this gives rise to the unnecessary additional constraint 
|A t /A s | > 1 for the existence of the non-Abelian phase. 
Since our method of calculating the index does not rely 
on this mapping, we find quite generally that the non- 
Abelian phase survives even if the gap closes at some 
value of |At/A s | and momentum |k|. This is illustrated 
in Fig. [l]and Fig. g] 

As is clear from Fig. [TJ the excitation gap, which is 
defined as the momentum-space global minimum of the 
lowest energy positive eigenvalue of the Hamiltonian in 
Eq. ([I]) , starts being non-zero for small values of the ratio 
A t /A s (we have chosen parameters such that the condi- 
tion V 2 > Aj + /i 2 is satisfied). The gap decreases with 
increasing values of A t /A s , eventually vanishing at some 
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FIG. 2: Plot shows how the ratio A t /A s at gap closing scales 
with the Zeeman splitting V z for a = 1, fi = 1, and A s = 0.5 
(in units of n — ft 2 /2m). Note that, since V z > /i 2 + A 2 for 
the range of parameters in the plot, the topological invariant 
Ci = 1 in spite of the gap closing. Inset shows the dependence 
on A s of the ratio (At/A 3 )o at which global minimum of the 
gap shifts from k x ,k y 7^ to the origin. For large A s , the 
global minimum is always at k x — k y = 0. Following values 
have been used: V z = 6.0, fi = 4.0, and a = 1.0 (in units of 
t] = ft 2 /2m) 

(parameter-dependent) value of this ratio. However, for 
larger value of A t /A s the gap opens again. One of the 
principal finding of this paper is that, even though the 
gap closes at some value of A t /A s (and |k| 7^ 0), the 
index C\ = 1 on both sides of the gap closing point. 
This indicates that the non-Abelian phase survives on 
both sides of the gap closing point on the A t /A s axis. 
Therefore, the closing of the gap with increasing values 
of A t is not associated with a topological quantum phase 
transition. In Fig. [lj for A t / A s larger than the value for 
which the gap closes, the gap also shows a discontinuity. 
This is because it discontinuously shifts from being at 
a nonzero k to k = 0. We find that, for a given set of 
values for a,V z ,fi, and A s , the gap lies at |k| 7^ for 
small values of the ratio A t /A s . Beyond a certain value 
of A t /A s (which is parameter-dependent and is typically 
larger than the value for which the gap closes), the gap 
shifts to |k| — and, as illustrated in Fig. [I] is discon- 
tinuous when this shift occurs. In Fig. [2j we show the 



various parameter dependencies of the two special values 
of At/ A s , the one at which the gap closes, and the one 
at which the gap shifts in the momentum space and is 
discontinuous. 

It is worth reiterating the surprising result stated 
above (which we confirm below also by explicit real space 
solutions of the BdG equations): even though the exci- 
tation gap of the superconducting system does vanish 
at a certain value of A t /A s and is non-zero on either 
side of this axis, the system has precisely the same non- 
Abelian topological properties on both sides of the gap 
closing point. This is a consequence of the fact that gap 
closing is a necessary but not a sufficient condition for 
having a topological phase transition. The supercon- 
ducting states on both sides of the gap closing point are 
characterized by order parameter defects carrying zero- 
energy non-Abelian Majorana fermion excitations. The 
only topological phase transition that takes place in this 
system happens when the chemical potential or the Zee- 
man splitting are tuned to values such that V z = /i 2 +A 2 . 

The integer computed in Eq. ^ is essentially the same 
as the TKNN number which gives rise to quantized Hall 
conductivity and topologically protected edge currents in 
quantum Hall systems. [25] There is a one-to-one corre- 
spondence between quantum Hall systems and supercon- 
ductors with broken time reversal symmetry - the BdG 
Hamiltonian characterizing the quasiparticles in the lat- 
ter also describes quantum Hall insulators, with the su- 
perconducting gap being replaced by the bulk insulating 
gap. We shall explicitly show the connection between C\ 
and the TKNN invariant in section IV. 



II. EXISTENCE OF MAJORANA MODES IN 
NCS WITH ZEEMAN SPLITTING 

In this section we shall explicitly derive an existence 
condition for the non-degenerate zero-energy Majorana 
modes in the vortex core of a NCS with Zeeman split- 
ting. In the absence of a Zeeman splitting, time reversal 
invariance guarantees two degenerate Majorana modes 
bound to vortex cores. We shall show that the existence 
criterion for single Majorana fermions bound to vortex 
cores is the same condition derived in the previous sec- 
tion. The BdG Hamiltonian in the presence of a vortex 
with winding number n can be written in real space using 
polar co-ordinates as 



H 



BdG 



[H t z + A s (r)r+ + h.c.]S(r - r') + - [A t (r,r')a + + A*(r,r')cx_] [cos(n6)T x + sm(n0)T y 



(4) 



In the above equation H = — h 2 ~S/ 2 /2m — /i + V z a z + tonian in the momentum space where the p-wave gap 
a(a x p).z. For compactness we shall write the Hamil- 
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is given by — iA t (a x p).z[cos(n8)T x + sm(n8)T y ]. The 
Hamiltonian Q can be diagonalized by noting that 
one can construct a pseudo-angular momentum operator 
J z = L z + (a z — nr z )/2 that commutes with the Hamilto- 
nian. One can therefore use a canonical transformation 
i?BdG = e -1 -'* H^dG^ to make the Hamiltonian inde- 
pendent of 9, where j z — rrij — u z /2 + nr z /2 and uij are 
the eigenvalues of the operator J z . We next note that 
under a particle- hole transformation mj — > —mj. Since 
we are interested in the Majorana solution we therefore 
set mj = 0. To further simplify the analysis we focus 
on vortices with unit winding number n — 1 and make 
another canonical transformation e^^^ifBdGe 10 * 7 ^ 4 
that makes the Hamiltonian real. The E = solu- 



where r] — ti 2 /(2m). 

Following p3] we next approximate the radial depen- 
dences of A s (r) and A t (r) by A Sjt (r) = for r < R, 
where R is the size of the vortex core and A Sj t(r) = A° t 
for r > R. Since both the order parameters are zero in- 
side the vortex core, one can construct explicit analytical 
solutions in terms of Bessel functions. The characteris- 
tic equation leads to two linearly independent solutions 
that correspond to the two Fermi surfaces obtained by 
the intersection of the two bands with the Fermi level. 
Out side the vortex, we seek solutions of the type: 

/ ut(r) \ _ e- / Pf (l/r) \ 

{ H(r) ) ~ r i/2 p x (l/r) J W 

such that Prj(x) are analytic functions of x. This leads to 
a convergent power-series solution in 1/r, whose zeroth 
order term can be obtained by setting 1/r = which 
leads to the following equation: 

f C + V, + zA t AA S -za \( Pt(0) \ n m 
\ -XA s + za £, z -V z +zA t ) \ Pi (Q)J ~ [ > 

where £ z = —r/z 2 — p. The two families of solutions 
for A = ±1 are related to each other by z — » — z. The 
values of z that are consistent with the above equation 
are obtained by solving 

(& ± zA t f - V z 2 + (za T A s ) 2 = (8) 

The sign of the product of the roots of this quartic equa- 
tion is given by S = sgn(n„(z„)) = sgn(A^ + p 2 — V 2 ). 
When A 2 + /i 2 — V 2 > there are 2 real roots and 2 
complex roots, whereas when A 2 + p? — V 2 < there are 
4 complex roots. The different possibilities arc enlisted 
in Table- 1. 

Figure (|3| show the roots of Eq. ^ for sgn(Co) = ± and 



tion of the resulting BdG Hamiltonian would in gen- 
eral yield two solutions that are complex conjugate pairs: 
^(r) and \?*(r). Also, due to the particle-hole symme- 
try of the BdG Hamiltonian, if W(r) is a solution then 
so is a y Ty^(r). If these two solutions are related by 
^(r) = \a y T y ty(r), then (cr y T y ) 2 — 1 implying A = ±1. 
As a final step, we use the relation u y T y = A to decou- 
ple the particle and hole sectors in the pairing terms by 
writing t x — i\a y r z . The pairing terms are then given by 
(i\a y T z )[A s — iA t (a x p).z] and the particle-hole sectors 
given by r z = ±1 are thus completely decoupled. 

The BdG equations for the n — 1 vortex with zero 
energy written in real space is then given by 



(5) 
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4 Complex roots 
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4 Complex roots 


For all 4 roots 




+1 


2 Real, 2 Complex roots 


For 1 root 
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2 Real, 2 Complex roots 


For 3 roots 



TABLE I: Summary of solutions for Eq. (|8| for sgn(Co) = ± 
and A = ±1. 



A = ±1 and a given set of parameters. Below we show 
that in the first three cases, there is no unique solution for 
^(r), whereas in the fourth case an unique zero-energy 
Majorana solution exists. 

In order for the solution ^ to be normalizable we need 
Re(z n ) < for r > R. Further, an unique solution for 
the wavefunction ^(r) of the zero-energy non-degenerate 
Majorana mode, would require the two-component wave- 
functions \&(r < R) and 4"(r > R) to satisfy four 
boundary conditions for \ff(r) and ^'(r) at the bound- 
ary r — R. One additional condition comes from the 
normalization condition for \I/(r). This leads to a to- 
tal of five equations for \l/(r) and hence an unique so- 
lution (|6| exists when there are exactly five unknown 
co-efficients in W(r) — ^ n c n ^ n (r) corresponding to the 
roots z n of Eq. Q. From the discussion in the previ- 
ous paragraph, we find that when A 2 + p 2 — V 2 < 
and A = +1, then *(r < R) = J2 n =i,2 c n^n(r) and 



-V (d 2 r + \d r ) +V Z — fi 
— XA s (r) — ad r 



A t d r XA s (r) + a (d r + ±0 r ) 

-n (d 2 + \d r - - v z - [I - A t (d r 



*„(r) = 



() 



a 3 - 

N 

1- 


* 
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— i 1 1 — e- 

-0.5 

-1- 
-3- 


1 1 1 — 

7 0.5 

- R * 

* 




c 3 - 

2- 

N ~i- 


* 


1 1 » ' 

-2 -1 

-1- 
-2- 
-3- 


' 1 1 1 * 

7 1 2 

* 



3- 

b _2- 




1 1 1 — 8- 

-0.5 
• "I" 
_2- 

* -3- 


1 1 1 

7 0.5 




d * 3 - 

N 1- 

» 1 1 >-0- 




-2 -1 

-1- 

-2- 
* 

-3- 


1 1 

7 1 2 



/j,(y > 0) = 0, and look for solutions of the type ^(y) = 
a n e~ ZnV . The zero energy mode of the Hamiltonian 
written in the reduced 2x2 space (say A = +1 channel) 
then takes the form 



Cm +V z + zA t XA S - z n a 
-AA S + z n a £ Zn -V z + z n A t 



u» = (10) 



where £ Zn = — rjz 2 — /i. This is the same equation as 
obtained for the solution outside the vortex core. Match- 
ing the boundary conditions and normalizing the solution 
yields five equations which bear an unique solution when 
the number of unknowns is also five. Following the anal- 
ysis in the previous section one can see that the latter 
condition is satisfied when V 2 > fj, 2 + A 2 . 



FIG. 3: Nature of roots for Eq. <|8J> for the following parame- 
ters: a = 1, A s = 1, A t = 1, and (a) fj, = 4, V z = 2, A = +1, 
(b) p = 4, V, = 2, X = -1, (c) fj, = 2, V z = 4, A = +1, (d) 
fj, = 2, V z = 4, A = — 1 (in units of 77 = fi 2 /2m) 

ty(r > R) = J2 n =i 2 3 c «^n( r ) an d hence there are 
exactly five unknown co- efficients. Hence, there is an 
unique zero-energy Majorana solution in this case. On 
the other hand, when A^ + fi 2 — V 2 > 0, one can easily 
check that either there are no normalizable solutions or 
the number of unknown co-efficients do not equal to the 
number of equations and hence no unique zero-energy 
solution exists. 



III. NON-CHIRAL EDGE MODES IN NCS 
WITH ZEEMAN SPLITTING 

It is well known that topological phases are charac- 
terized by the existence of gapless edge states. Here we 
shall explicitly demonstrate the existence of edge states 
in NCS with Zeeman splitting and how the edge states 
are related to the existence of Majorana modes in vortex 
cores discussed in the previous section. Without loss of 
generality we consider an edge perpendicular to y. The 
Hamiltonian, which is now a function of k x can be shown 
to be symmetric under k x — > — k x . This property of the 
Hamiltonian follows from particle-hole symmetry. There- 
fore, in analogy with the calculation of the Majorana 
modes we set k x — 0. The resulting Hamiltonian is then 
a single band ID BdG Hamiltonian with spin-orbit cou- 
pling given by 

H B dG = [-rjdy - fi(y) + V z a z - iaa x d y ]r z 

+(A s {ij)t + + h.c.) + A t (y)(-ia x d y )(ia y T z ) (9) 

We note that the non-degenerate spinor solution has the 
Majorana form , F = (u, ia y u*), where u is a two spinor. 
Writing u = ur + iui, one obtains a pair of equations 
in which ur and ui are decoupled and are completely 
analogous to the A = ±1 channels for the Majorana case. 
One can define the edge at y = as ji(y < 0) = /j, and 



IV. PFAFFIAN Z 2 INVARIANT AND CHERN 
NUMBER INVARIANT FOR BDG 
H AMILTONI ANS . 

In the previous sections, we characterized the time- 
reversal symmetry breaking non-centrosymmetric super- 
conductors first by the first Chern number topological in- 
variant Ci and then by the sign P — sgn(/i 2 + A 2 — Vf ). 
The first Chern number C\ is an integer (Z) invariant 
while the latter quantity is a invariant that charac- 
terizes whether the system is in a phase that supports 
Majorana fermions or not at defects. Numerical calcu- 
lation of C\ suggests that both P = ±1 and C\ = 0, 1 
are consistent characterization of when the system sup- 
ports a Majorana Fcrmion. In this section, we explicitly 
show that P is a special case of a Pfaffian invariant for 
Hsdcik = 0), which is shown to be equal to the parity 
of the first Chern number C%, analogous to the case of 
topological insulators [2T)] . 

A. Pfaffian Invariant 

The k = (k x ,k y ) dependent Bogoliubov-de Gennes 
Hamiltonian Hb,ig(^) has a well-known particle-hole 
symmetry which is written as 

EH BdG (k)E- x = AH* BdG (k)A = -H BdG {-k) (11) 

where 5 = AK and A = (<J y T y ) and K is the complex 
conjugation operator. Here o x ,%t,z are the Pauli spin ma- 
trices and T, j , are the Nambu particle-hole matrices. 

The particle-hole symmetry S, similar to the time- 
reversal symmetry, maps fe — » — fc. Therefore similar to 
time-reversal invariant topological insulators, some of the 
topological properties of particle-hole symmetric BdG 
Hamiltonians can be extracted from the Hamiltonian at 
particle-hole symmetric k-points such that k = —k + G 
where G is a reciprocal lattice vector. We refer to such 
points as K. The BdG Hamiltonian at such points 
Hb<ig(K) is explicitly particle-hole symmetric. Given 
such a Hamiltonian, we can define a matrix W(K) which 
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can be shown to be anti-symmetric by virtue of the 
particle-hole symmetry as follows: 

W T {K) = A T H T (K) = —AAH(K)A = -W(K). (12) 

Moreover one can show that i n Pf(W) is real since 

Pf(W(K)) = Pf(W(K))' = Pf(H" Bia (K)A) 

= Pf(W T (K)) = (-l) n Pf(W(K)) 

where H BdG {K) is a 2n x 2n Hermitian matrix. Using 
this one can define the sign 

Q(H BdG (K)) = S gn(i n Pf(H BdG (K)A)) (13) 

for any particle-hole symmetric BdG Hamiltonian. 
The function Q defined above can change only if 
P f (H BdG {K)A) vanishes. Since 

Pf(H BdG (K)A) 2 = Det(H BdG {K)A) = Det(H BdG ), 

this implies that Q can only change sign where the gap 
of the BdG Hamiltonian vanishes at K. Therefore Q 
defines a topological invariant for the space of particle- 
hole symmetric BdG Hamiltonians. Note also that this 
topological invariant does not rely on any other sym- 
metry such as time-reversal symmetry. The invariant Q 
is related to the invariant suggested by Kitaev for one- 
dimensional p x superconductors [5] . It can be easily ver- 
ified that for the Hamiltonian in Eq.Q, the invariant 
Q(H(0))= S gmTj-A2- M 2). 

The Pfaffian invariant Q can be computed in terms 
of the eigenvectors of H BdG (K). Suppose H BdG (K) is 
diagonalized by the transformation 

H BdG (K) = U(K)D(K)U\K) 

where D(K) is a diagonal matrix of eigenvalues ordered 
in descending order of value. The columns of the unitary 
matrix U(K) are the eigenvectors of H BdG (K). The 
positive energy eigenvectors in U are chosen to be re- 
lated to the negative energy eigenvectors by particle-hole 
symmetry. With such a convention for the eigenvectors 
and eigenvalues, the unitary matrix U(K) satisfies the 
constraint 



AU(K) = U*(K)T 



where 



r = 



In 
In 



and /„ is the n x n identity matrix. The Pfafnan of 
H BdG {K)A is then 

Pf(H BdG (K)A) = Det(U(K))Pf(D(K)T) (14) 



where Pf(D(K)T) = ]J n>0 E n (K). Therefore the 
topological invariant Q is determined by the sign of 
Det(U(K)) as 



Q(H BdG {K)) = Det(U(K)). 



(15) 



The above formula can be used to evaluate 
Q{PiBdG{K)) at the particle- hole invariant points 
K = (b,±f),(±f,0),(±f,±S) using the fact 
that tK 3> Vz,A s ,/i. To dominant order in 
a, Hb<ig{K) — ejcTz- Therefore at these points 
Q(H BdG (K)) = 1. 



B. Relation to Berry connection. 

The evolution of the eigenvector matrix U (k x , k y ) with 
k x can be computed by introducing the multi-band vector 
potential 

A m . n {k x , ky) = i(u m {kxi ky)\V\u n (k x , ky)) . 

The evolution of the eigenvector matrix can be written 
in terms of this potential as 

d kx U(k x ,k y )=iU(k x )A^(k x ,k y ) 

U (k x , ky r^kUik^ky) = iA^(k x ,ky) 
d k JOg{U(k x ,ky))=iA^\k x ,k y ) 

d ksc \og(Det(U{k x , k y ))) = iTr(A (x \k x , k v )) 



P{ky = 0,tt) 



Det(U(k x = 7T . ky ) ) _ inS 



Det(U(k x = 0, k y )) 
Here S is integral of the trace and is calculated as 

S= - [ dk x Tr(A^(k x ,k y )) 
^ Jo 

r * dk. 



This integral by itself does not reduce to a BZ integral 
since the integral over k x extends only to 7r. Using parti- 
cle hole symmetry of the BdG Hamiltonian it is possible 
to pick Bloch states at — k x such that 

u- m (-k x ,k y =0,7r) = a y T y u* m (k x ,ky = 0,7r) (16) 

where —to is the state in the band —to. In this gauge, 
the connection away from k x = ir must satisfy 

^mikx^ky) d kx u m {k x ^ ky) u, m (k x: ky) d kx u m (k x: ky) 

= ~ M — m( — k x , ky)df. x U— m {— k x , ky) 

= "-mC*!! ky)d k >^U— m (k x , ky) \k' x =—k x ■ 



Therefore the integral of the trace S becomes 
S = I dk x S (u m (k x , ky ) \dk x \u m (fe X) ky )) 

■ / ° m<0 

^ ] ( M m (fexi fey ) l^fc x (fea ; ky)} 
m>0 

r 

d^x ^ ^ ( u m(k X i k y )\dk :l .\u m (k x , ky)) 

m<0 

dk x (u m (k x , ky)\df t:r \u rn {k x , ky)) 

T m<0 
: 

dk S ' (u m (fej; , fcy) l^fe^ |w m (k x , ky)) . 



m<0 

The integral S can be written in terms of a t/(l) vector 
potential a(k x ,k y ) 

— -a x (k x ,k y ) (17) 



where 



a(k x ,k y ) = }^ (u m (k x , ky)\iV\u m (k x , ky)). (18) 
m<0 

The line integral over fc x can be converted to a loop 
integral that encloses a finite area of the BZ torus by 
subtracting the integrals at k y = and k y — tt i.e. 

5 - S(k y = tt) - S(jfc„ = 0) 

x&xij^xi ky — tt) I dk x o, x [k x •> ky — 0) 



dk x Ojx{kxi ky 



(19) 



Using Stokes theorem this integral can be written in 
terms of the Berry curvature f xy = df. x a,k y — dk y 0Lk x as 



dk x a X (k X ,ky) - j dk X dkyf X y{k Xl ky) (20) 

which as pointed out by TKNN [28l [30] is written as 
dk x / dky y)(fm ~ fn) 

-TT JO m „ 



(U m 1 <9fc x H B dG I Un) (u n \ 8 ky H B dG I 

(E m — E n ) 2 



(21) 



where f m ,n are the Fermi occupation functions. 

Using particle-hole symmetry, the above equation can 
be extended to the entire BZ to yield the final result 



P 



where S = j BZ d 2 kf xy and 

\d kn .H B dG\u n ) (u n \dk y HBdG\ 
fxy = / \fm ~ Jn)~ TTp _ Tp \2 ~~ ' 

(23) 

Since, as argued in the previous sub-section, the Q 
invariant at the particle-hole symmetric points other than 
K = are 1, the results of this sub-section shows that the 
Pfaffian topological invariant Q(H B dc{k x =0,k y = 0)) 
is related to the TKNN number S by 

Q(H BdG (k x = 0, k y = 0)) = e" s . (24) 

This is the central result of this section. 

C. Connection to Green function form of invariant. 

The TKNN invariant can be expressed compactly in 
terms of the Green function as [27] 

Ci = J d 2 kdu}Tr[Gd kj; G- 1 Gd ky G- 1 Gd UJ G- 1 } 

- Tr[Gd ky G~ 1 Gd k]: G~ 1 Gd UJ G~ 1 } (25) 

where G^ 1 (k x , k y ,uS) = cj — H(k x , ky). To see that this is 
indeed the same as the TKNN expression, we note first 
that the above equation is simply 

Ci = J d 2 kdujTr[G 2 d k!c G- 1 Gd kv G- 1 ] 

-Tr^dk^Gd^G- 1 }. (26) 

Expanding the Green function in terms of eigenstates we 
see that the frequency integral in the momentum can be 
written as 



duji y] 7 ^r^iumdk^Hlun) — (u n \d k H\u Tl 

±-s (w - E m y u - E n 



[fm fn) 



(u m dk x H\u n }(u n \dk v H\u m ) = fxy. 



Q(H B dG{k x = 0, ky = 0))Q(H BdG (k x = 7r, k y = tt)) 
~~ Q{H BdG {k x = n, k y = 0))Q(H BdG (k x =0,k y = tt)) 
= e" s (22) 



The invariant Ci then becomes the TKNN invariant 
C t = S= f «/-/,•/, „■:/.•,./.•,). (27) 

J BZ 

V. CONCLUSION 

We have considered non-centrosymmetric supercon- 
ductors in the presence of a Zeeman splitting normal to 
the surface. The resulting two-dimensional superconduc- 
tor on the surface has both s-wave and p-wave pairing. 
We have derived the condition for which the system is in 
a non-Abelian phase by constructing a topological invari- 
ant. We find that the condition for a non-Abelian phase 
to exist is completely independent of the triplet pairing 
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potential. This is in contrast to the case of time reversal 
symmetric NCS, where a non-trivial topological phase 
exists only when the triplet pairing is larger than the 
singlet one [3T] . We find that even though the excitation 
gap of the superconducting system vanishes at a certain 
value of At/ A s and is non-zero on either side of this axis, 
the system has precisely the same non-Abelian topologi- 
cal properties on both sides of the gap closing point. The 
superconducting states on both sides of the gap closing 
point are characterized by order parameter defects car- 
rying zero-energy non-Abelian Majorana fermion excita- 
tions. There is no quantum phase transition, not even 
of the topological type, even though the excitation gap 
vanishes as the parameter A t /A s is tuned! The only 
topological phase transition that takes place in this sys- 
tem happens when the chemical potential or the Zeeman 
splitting are tuned to values such that V 2 = fi 2 + A 2 . For 
V z smaller than the critical value V£ = ^ fj, 2 + A 2 , the 
system is a regular non-centrosymmetric superconductor 
with a mixture of s-wave and p-wave pairing amplitude. 
On the other hand, for V z greater than the critical value 
V£, the system is a non-Abelian non-centrosymmetric su- 
perconductor with order parameter defects such as vor- 
tices carrying non-degenerate Majorana modes. 

We next look for non-degenerate Majorana solutions in 
vortex cores in non-centrosymmetric superconductors by 
explicitly solving the relevant BdG equations. We find 
that the same condition on the Zeeman splitting, chem- 
ical potential, and (only) the s-wave component of the 
order parameter needs to be satisfied for a unique Majo- 
rana solution to exist. As a corollary of this calculation, 
we find that the existence of the zero-energy Majorana 
modes in the bulk defects such as vortices is intricately 
related to the existence of gapless edge modes. The ex- 
istence condition for the non-Abelian state involves the 
quantity (V 2 — /i 2 — A 2 ), which is really the Pfaffian of 
the BdG Hamiltonian at momentum |k| = 0. Since the 
p-wave pairing amplitude in the Hamiltonian does not 
contribute to this Pfaffian, we arrive at the important 
result that the non-centrosymmetric superconductors are 



always in the non-Abelian phase as long as A s satisfies 
V 2 > [i 2 + A 2 , irrespective of the value of the p-wave 
component of the order parameter. 

We have derived the important result that the non- 
Abelian properties of the non-centrosymmetric supercon- 
ductors are completely insensitive to the magnitude of 
the p-wave piece of the order parameter by seemingly 
two independent methods: by using a Z topological in- 
variant and by the direct solutions of the BdG equations 
which gives us a condition involving the Pfaffian of the 
BdG Hamiltonian at k — 0. In the last part of the paper, 
we make the connection between the topological invari- 
ant, which is nothing but the TKNN number, and the 
Pfaffian of the ID BdG Hamiltonian for the Majorana 
mode. The Pfaffian is shown to be related to the first 
Chern number through the Berry connection. The fact 
that the superconductivity is intrinsic {i.e., it does not 
need to be induced by the proximity effect) and the non- 
Abelian phase can exist for any value of the triplet pairing 
amplitude makes non-centrosymmetric superconductors 
a promising candidate for the realization of Majorana 
bound states and non-Abelian statistics. 
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